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An unexpected insensitivity of the Fermi surface to impurity scattering is found in Ru substi- 
tuted BaFe2As2 from first-principles theory, offering a natural explanation of the unusual resilience 
of transport and superconductivity to a high level of disordered substitution in this material. This 
robustness is shown to originate from a coherent interference of correlated on-site and inter-site im- 
purity scattering, similar in spirit to the microscopic mechanism of superdiffusion in one dimension. 
Our result also demonstrates a strong substitution dependence of the Fermi surface and carrier 
concentration, and provides a natural resolution to current discrepancies in recent photoelectron 
spectroscopy. These effects offer a natural explanation of the diminishing long-range magnetic, 
orbital and superconducting order with high substitution. 



PACS numbers: 74.70.-b, 71.15.-m, 71.18.+y, 71.23.-k 

Chemical substitution is the most widely employed 
technique to induce high-temperature superconductiv- 
ity pj[5]. The most obvious effects of substitution are 
carrier doping and chemical pressure. On the other hand, 
such substitution should also introduce disordered impu- 
rity scattering, an effect mostly unexplored in the studies 
of electronic structure of real materials due to its com- 
plexity. Recently, owing to an advances in first-principles 
theoretical methods [3j [6], various novel and surpris- 
ing physical effects of disorder were found in the new 
Fe-based superconductors containing transition metal 
dopants [8-10 and Fe vacancies [11 . Naturally, one 
would wonder whether there are also unexpected disorder 
effects in isovalent substitution. An ideal candidate to in- 
vestigate is the Ru substituted Fe-based superconductor 
Ba(Fei_ a ,Ru :c )2As2(Rul22) [T2H2Q]. Unlike other exist- 
ing cases with isovalent substitution at the anion sites, 
Ru atoms replaces the most essential Fe atoms. Given 
that Ru is a 4d element, it must be quite different chem- 
ically from 3d Fe and therefore must introduce a strong 
impurity scattering potential. 

One obvious physical puzzle regarding Rul22 is the 
resilience of its superconductivity and transport against 
large concentration of disordered impurities. Indeed, 
samples are found to retain their superconductivity even 
with 40% Ru substitution of Fe [13] [15] . Given the cur- 
rent proposal that the superconducting order parameter 
is most likely of s± symmetry [4] , it is extremely puzzling 
how a sign-changing order parameter can survive such a 
large concentration of strong impurities. Similarly, given 
that Ru is more distinct from Fe than Co is, it is quite un- 
expected that Ru substituted samples exhibit a residual 
resistivity comparable to 8% Co substituted system [5] at 
a much higher 35% substitution level. 

Another important current issue is the Ru substitu- 



tion dependence of the electronic structure. Qualitatively 
similar to the Co substituted systems, Ru substitution at 
high enough level systematically suppresses the magnetic, 
orbital order, and then superconductivity. For Co substi- 
tuted systems, this behavior can be understood from the 
weakening of the nesting between the electron and hole 
pockets, since the former grow in size while the latter 
shrink, responding to the additional doped carriers [TP] . 
For Ru substitution, however, the size of the electron and 
hole pockets should remain balanced (and approximately 
nested), since Ru is supposed to be isovalent to Fe. It is 
thus not as obvious what substitution-dependent features 
in the electronic structure suppress the long-range order 
in this case. 

Currently, this issue is quite controversial in the 
field, due to seemingly contradictory experimental ob- 
servations. Some angular resolved photoemission spec- 
troscopy (ARPES) experiments [15] reported a nearly 
substitution-independent Fermi surface, and suggested 
that superconductivity emerges from dilution of the mag- 
netism, in contrast to other doped 122 systems. Other 
ARPES measurements, on the other hand, reported a 
large increase of the number of carriers [14 with in- 
creasing Ru concentration and a crossover from two- 
dimensional to three-dimensional structure of some of the 
hole-like Fermi surfaces [16]. Similar contradictions also 
emerge in current first-principles computations. Some 
concluded that Ru substitution induces no doping [18] 
and no changes in the low energy band structure [19], 
while others [20] found large changes. It is thus timely to 
investigate the substitution-dependence of the electronic 
structure of Rul22 with a proper account of the disorder 
effects to resolve the current debates, and if possible to of- 
fer an explanation for the suppression of long-range mag- 
netic, orbital and superconducting order upon increased 
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substitution levels. 

In this Letter, we address these important issues by 
studying the electronic structure of Ba(Fei_ cc Ru :c )2As2 
over the full range of substitution (0 < x < 100) tak- 
ing into account the realistic disorder effects in first- 
principles calculation. This is made possible via the re- 
cently developed Wannier function based effective Hamil- 
tonian method for disordered systems [6 . Surprisingly, 
while large scattering is found on the entire Fe band com- 
plex, for all substitution levels, the states near the chem- 
ical potential remain very sharp and coherent. This un- 
expected insensitivity of the quasiparticles to impurity 
scattering is here traced back to a coherent interference 
between on-site and off-site impurity effects. In spirit, 
this is very similar to the microscopic mechanism that 
gives rise to superdiffusion in ID theory [21J, and can 
thus be considered as the first realization of this exotic 
phenomenon in real materials. 

This insensitivity to impurity scattering provides a 
natural explanation of the amazing resilience of transport 
and superconductivity in the presence of a high level of Fe 
site substitution in this specific material. In addition, our 
results reproduce the measured spectral functions and 
resolve the current controversy in the interpretations of 
the Ru substitution dependence. We find a systematic 
reduction of carrier density upon Ru substitution accom- 
panied by an enhanced 3D character of the Fermi surface. 
These effects lead naturally to suppression of long-range 
magnetic and superconducting order at high substitution 
level. These findings highlight the general need to in- 
corporate disorder effects in most chemically substituted 
systems. 

We first evaluate the configuration-averaged spectral 
function (A n (k,uj)} of Wannier orbital n, frequency co>, 
and crystal momentum within the approach of Ref. 
El and with the help of the band structure unfolding 
method [3J. The resulting (A n (k,uj)} is found con- 
verged within the resolution of the plots after averag- 
ing over 10 large randomly-shaped supercells contain- 
ing 400 atoms on average with random Ru substitu- 
tions. The influence of Ru substitution on the Hamil- 
tonian is extracted from the results of several den- 
sity function theory (DFT) calculations: the undoped 
BaFe2As2 or BaRu2As2, as well as the impurity super- 
cells of Ba2FesRuAs4 or Ba2FeRusAs4. The low energy 
Hilbert space is taken within [-6,4] eV consisting of Wan- 
nier orbitals of Fe-d, Ru-d and As-p characters. Careful 
benchmarks of the quality of our effective Hamiltonian 
are conducted against full DFT calculations using stoi- 
chiometric ordered "impurity lattice" test cases [22] . The 
entire study is conducted using the 2-Fe Brillouin zone 
required by the symmetry [23 , and unfolded to the 1-Fe 
zone [3, 23] when direct comparison with ARPES spectra 
is needed. 

Figure [4] summarizes our resulting spectral func- 
tions over the entire range of substitution for 



Ba(Fei_ a ,Ru :c )2As2 (x = 0.0,0.21,0.38,0.55,0.75 and 
1.0). A clear broadening of quasiparticle spectral lines 
in both momentum and frequency is observed over the 
entire Fe d-band complex, reflecting the finite mean free 
path and lifetime of carriers in these states due to dis- 
order. The broadening is naturally the strongest for the 
x = 0.55 case, since it is the most disordered one shown. 
(As x approaches or 1, the system starts to become 
cleaner, approaching the pure end compounds.) Such a 
strong impurity scattering is expected, since the 4d level 
of Ru is about 0.8 eV lower in energy than the 3d level 
of Fe, and the different spatial extent of the orbitals also 
implies large changes in their hopping strength. 

Strikingly, the bands near the Fermi level (shown in the 
second panel of Fig. 1) remain very sharp, as if carriers 
occupying states near the Fermi level don't scatter from 
the impurity. Correspondingly, the configurational aver- 
aged spectral function [Fig. 2(a)] exhibits a very sharp 
peak at the Fermi level with a well-defined dispersion. 
(In comparison, around ~ — 1 eV the spectral function 
shows a much broader peak, reflecting strong impurity 
scattering effects.) Therefore, the puzzling resilience of 
transport and superconductivity against high substitu- 
tion level in Rul22 is now understandable simply from 
the lack of net impurity scattering for states near the 
Fermi level. 

Fig. [5jb)-(d) give further insights into the microscopic 
origin of this unusual insensitivity to disorder near the 
Fermi level. Here we compare spectral functions under 
the influence of different partial components of the im- 
purity potential in the Wannier function basis. With 
only the diagonal impurity potential (on-site disorder), 
Fig. [5^c) shows a strong smearing of the entire d-band 
complex, including states near the Fermi level. On the 
other hand, with only the off-diagonal impurity poten- 
tial (off-site disorder), Fig. [5jd) shows a well structured 
smearing of the band structure: the farther away from 
the center of the d-bands, the more the states are af- 
fected by the impurity. In particular, near the center 
of the band, the effects of the off-site impurity poten- 
tial diminishes. This is understandable since the off-site 
terms of the Hamiltonian are responsible for the band 
dispersion, and their effects is maximum at the band 
edges. One would expect a similar structure for the effec- 
tiveness of off-site impurity potential, particularly posi- 
tive/negative above/below the center of the band. (The 
Ru orbital is bigger than Fe, so the impurity potential 
tends to enhance the band width.) A similar conclusion 
can be reached by imagining the band structure interpo- 
lated between pure BaFe2As2 and pure BaRu2As2- 

Upon combining both on-site and off-site disorder po- 
tential, Fig. [5jb) shows that the energy range with weak 
impurity scattering moves up, close to the Fermi level. 
This can now be understood as a consequence of the 
cancelation of (negative) on-site and (positive) off-site 
impurity potential above the center of the band. Indeed, 
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(a)x=0 (b)x=0.21 (c)x=0.38 (d)x=0.55 (e)x=0.75 (f)x=1.0 

a=3.96 c=13.01 z =0.35 a=4.00 c=12.85 z =0.35 a=4.03 c=12.75 z =0.35 a=4.06 c=12.59 z =0.35 a=4.10 c=12.44 z =0.35 a=4.15 c=12.25 z a =0.35 
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FIG. 1: (color online) The lattice constants (A) [4j [6j [13], band structure (the first panel), low energy band structure(the 
second panel) , the Fermi surface around Y point with k z = (the third panel) and around Z point with k z = tt (the last panel) 
for Ba(Fei_ x Ru x ) 2 As 2 with different Ru concentration x = (a), 0.21 (b). 0.38 (c), 0.55 (d), 0.75 (e) and 1.0 (f). Orbital 
characters of Fe, Ru, and As are distinguished by color given in the legend. Resulting carrier density per Fe, and density of 
states at the Fermi level N(Ef) per (eV Fe) are given below the plots. 



Fig. [5|e) gives an example of the net coupling involving 
8th/7th and 9th/10th bands near the T point resulting 
from a single impurity for all three cases. One sees that 
near the Fermi level, the impurity scattering between 9th 
and 10th bands has large on-site contribution (similar to 
that between 8th and 7th bands), but it is almost en- 
tirely canceled by the equally large off-site scattering. It 
is important to note that such coherent interference can 
only take place because the positions of on-site and off- 
site impurity potential are always correlated: both are 
associated with the location of Ru substitution. There- 
fore, moving the impurity position only adds the same 
overall phase to both of the terms but does not affect 
their interference. 

This effect is in spirit very similar to the microscopic 
mechanism that gives rise to the unusual superdiffusion 
in ID models [21 j. Due to the poor connection between 



sites, (non-interacting) ID systems are known to be ex- 
tremely sensitive to disorder. Any tiny on-site disorder 
leads immediately to a total localization of the system. 
Surprisingly, it can be shown [21] that the system can 
still overcome the localization to give a "superdiffusive" 
metallic transport if certain special correlations exists 
in the position of the impurities to ensure a coherent 
interference of their scattering. What we find here in 
Rul22 is probably a more general (and realistic) case, in 
which the interference takes place at a single impurity, 
but between on-site and off-site contributions. From this 
point of view, Rul22 can thus be considered the first re- 
alization of such an exotic superdiffusion mechanism in 
real materials. This realization manifests itself with a 
significant consequence that superconductivity survives 
to a high substitution level, due to the insensitivity of 
states near the chemical potential to the impurity scat- 



4 



(a) <A(r»> (eV 1 ) 



(b) x=0.38 



(c) diagonal disorder (d) off-diagonal disorder 



I 



(a) x=0.02 



(b) x=0.21 



(c) x=0.36 



(d) x=0.55 



(e) 



10 20 30 r M X TZRAZTMX TZRAZTMX TZRAZ 




k 


j 


k 


j 


<kj lAlkj > 


<kjlA D lkj> 


<kjlA 0D lkj> 


r 


9 


0.3rM 


10 


-0.0578 


-0.3023 


0.2445 


r 


8 


0.3rM 


7 


-0.2938 


-0.3316 


0.0378 



FIG. 2: (color online) Spectral function (a) and band struc- 
ture for system with two types of disorder (b), with only 
the diagonal disorder (c), and with only the off-diagonal 
disorder(d) for disordered Ba(Feo.62Ruo.3s)2As2. (e) single- 
impurity induced coupling between the Bloch orbitals near 
and far from the chemical potential, indicated by the black 
and red segments in (d). 



FIG. 3: (color online) Comparison of the Fermi surface at 
k z = tt plane between the ARPES results [15] (top panel) 
and our results (bottom panel) in 1-Fe zone for different Ru 
concentration x = 0.02 (a), 0.21 (b), 0.36 (c) and 0.55 (d). 
The dashed lines marks the true 2-Fe Brillouin zone bound- 



tering. In general, similar interference should occur in 
some other materials with disordered impurities as well, 
since realistic impurity induced on-site and off-site ef- 
fects are always correlated in position. This might also 
explain why metallic transport can be found in some thin 
nanowires [26] [27] and quasi- ID materials [28] [29] , de- 
spite the unavoidable occurrence of disorder. 

Next, considering Ru substitution dependence of the 
electronic structure, Fig. 1 show several clear trends. 
First, the hole pockets in the center of the Brillouin zone 
becomes more three-dimensional, shrinking significantly 
in size in the k z = plane and eventually disappearing, 
while surviving in the k z = tt plane. Specifically, the 
pocket with large z 2 character shows the strongest 3D 
dispersion, growing in size in the k z — tt plane. On the 
other hand, the electron pockets near the corner of the 
zone shows only a slight substitution dependence. This 
large change is quite consistent with the recent ARPES 
observation [16] , but seems to contradict the earlier ob- 
servation of substitution independent Fermi surface [T5] . 

This discrepancy can, in fact, be easily resolved by 
accounting for the matrix element of the incident photon 
in the experiment [30 . The photon polarization used in 
Ref. [15] is along the y direction, which couples mostly 
to the yz orbitals. Figure 3 shows the yz character of 
our theoretical unfolded band structure [3] , showing little 
change upon Ru substitution at the k z = tt plane, in 
excellent agreement with the experiment. 

The second physical effect of Ru substitution revealed 
in Fig. 1 is an physical reduction of carrier density, a 
real doping effect, and correspondingly a reduction of 
the density of states (DOS) at the chemical potential. 
In typical systems with only one sign of carriers, carrier 
doping is typically achieved via substitution of different 
valence. Here, Ru substitution, while typically consid- 
ered isovalent, actually reduces the carrier concentration 



significantly, from 0.16/Fe at x = to 0.12/Fe at x = 1, 
about a 25% reduction (c.f. bottom of Fig. 1.) This 
takes places through reducing the 3D volume of both the 
electron and the hole pockets (while keeping charge neu- 
trality), as if the system is being drained of both electron 
and hole simultaneously. Microscopically, this originates 
from a larger splitting between the conduction and va- 
lence bands, due to the enhanced hopping via Ru sub- 
stituted sites. Such chemical "pressure" effects should 
survive the mass renormalization observed in real mate- 
rials. 

Our result, which is perfectly consistent with the recent 
ARPES measurement [16] , suggests a more complex pic- 
ture. Indeed, we found that the green z 2 pocket increases 
in size at kz = tt upon Ru substitution, as observed ex- 
perimentally [16]. However, the total volume of the hole 
pockets (or equivalently that of the electron pockets) ac- 
tually reduces, as reflected in the corresponding carrier 
density shown at the bottom of Fig. 1. 

The above substitution dependence of the electronic 
structure corresponds nicely to the overall features in the 
phase diagram of this material, from the standard weak 
coupling perspective. Both the enhanced 3D nature and 
the reduced DOS tend to suppress the long-range mag- 
neic/orbital order, and eventually the superconductivity 
with sufficient high substitution level (the "overdoped" 
regime.) Indeed, the enhanced 3D band structure weak- 
ens the nesting condition for magnetic/orbital correla- 
tion, and presumably the glue for superconductivity in 
the "overdoped" regime. Similarly, the smaller DOS re- 
duces the Fermi surface instability in both the magnetic, 
orbital, and pairing channels. In essence, the case of Ru 
substitution offers a unique test case with strong substi- 
tution dependence of the electronic structure but with- 
out tipping the balance between electron and hole pock- 
ets. This unique feature of Rul22 should serve as a good 
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qualitative benchmark for all proposed theories of super- 
conductivity. 

In summary, a remarkable insensitivity of the states 
near the chemical potential to disordered impurity scat- 
tering is found in Ba(Fei_ a; Ru :c )2As2, despite the strong 
impurity potential of Ru. This offers a natural expla- 
nation of the unusual resilience of transport and super- 
conductivity against a high impurity level. This exotic 
behavior originates from a coherent interference of on- 
site and off-site disorder potential, and thus can be con- 
sidered the first realization of super-diffusion mechanism 
in real materials. Our result, while being very consis- 
tent with existing ARPES measurement, resolves a few 
current discrepancies in the field, and leads to new phys- 
ical conclusions. Specifically, Ru substitution is found 
to physically reduce the carrier density, instead of being 
"iso-valent" . Our findings offer a natural explanation of 
the diminishing long-range magnetic, orbital and super- 
conducting order at high substitution level. 
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of Energy, Office of Basic Energy Sciences DE-AC02- 
98CH10886, DE-FG02-05ER46236 and DOE-CMCSN. 
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I 

[SUPPLEMENTARY INFORMATION] 

Details of the density functional theory calculations 

We investigate the electronic structure of Ba(Fei_ ;E Ru :E )2As2 for the six doping levels listed in the table at the 
right bottom of figure |4j For the end points x = and x = 1, ordinary Density Functional Theory (DFT) is used, 
while for the doping levels in between the effective Hamiltonian method [2] with the help of the band structure 
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FIG. 4: (Left panel) Undoped normal cell BaFe2As2/BaRu2As2 , and (right panel) single Ru/Fe impurity supercell 
Ba2FesRuAs4/Ba2FeRu3As4, and (right bottom) lattice constants used for the different doping levels. 



unfolding method [3 is employed to be able to properly treat the influence of disorder. As an input for the effective 
Hamiltonian method [2] we perform two DFT calculations for each of the four intermediate doping levels. For x = 0.21 
and x = 0.38 we use BaFe2As2 as the undoped system and treat Ru as an impurity, while for x = 0.55 and x = 0.75 we 
use BaRu2As2 as the undoped system and treat Fe as an impurity. The lattice structures of the undoped normal cells, 
BaFe2As2/BaRu2As2 , are depicted in the left panel of figure [4] The primitive normal cell lattice vectors expressed 
in Cartesian coordinates are given by: a\ = — \ax-\-\ay-\-\cz ; <22 = \ax — \ay +\cz ; as = \ax + \ay — \cz . 
The lattice parameters for the doping levels x = 0,x = 0.38,x = l are taken from[3H6], respectively, and those 
for the doping levels x = 0.21, x = 0.55, x = 0.75 are obtained from linear interpolation between x = and 
x = 1. To capture the influence of the Ru/Fe impurity a two times larger supercell is used to properly treat the 
non-local influence on the nearest As and Fe/Ru orbitals. The corresponding Ba2FesRuAs4/Ba2FeRu3As4 supercells 
are depicted on the right panel of figure [4] and its primitive lattice vectors expressed in normal cell lattice vectors 
are given by: A\ = ci2 + as ; A2 = a\ + as ; As = a\ + a^. We applied the WIEN2K[1 implementation of the full 
potential linearized augmented plane wave method in the local density approximation. The k-point mesh was taken 
to be llxllxll for the undoped normal cell and 17x17x5 for the supercells respectively. The basis set sizes were 
determined by RKmax=7. 



Benchmarks of the effective Hamiltonian against DFT 



To explore the accuracy and efficiency of the effective Hamiltonian method [2 for the case of Ba(Fei_ :E Ru ;E )2As2, 
we present comparisons of spectral functions A n (k,uS) calculated from the full DFT and the effective Hamiltonian (see 
figures [5][6| . The size of the deviations between the full DFT and the effective Hamiltonian should be compared with 
the size of the impurity induced changes. For this purpose the spectral function of the undoped BaFe2Se2/BaFe2Ru2 
is also plotted as a reference for each benchmark. The Ba4Fe4Ru4As8 supercell in figures [5] and [6] is designed to test 
the linearity of the impurity influence (see formula (1) of Ref. [2]) and the partitioning of the impurity influence from 
its super images (see section IV of the supplementary material of Ref. [2 j). The basis set of Linear Augmented Plane 
Waves (LAPW's) used in the full DFT is ~ 30 times larger then the basis set of Wannier functions used in the effective 
Hamiltonian method. Since the number of floating point operations of diagonalization depends cubically on the size of 
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the matrix this implies an efficiency increase by a factor of 30 3 ~ 3- 10 4 . Furthermore the full DFT calculation involves 
multiple self consistent cycles (typically 15) whereas the effective Hamiltonian method (as currently implemented) 
requires only a single diagonalization, which increases the efficiency by another order of magnitude to ~ 6 • 10 5 . 
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